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Abstract 



This paper focuses on isochronicity of linear center perturbed by 
a polynomial. Isochronicity of a linear center perturbed by a degree 
four and degree five polynomials is studied, several new isochronous 
^ I centers are found. For homogeneous isochronous perturbations, a first 

integral and a linearizing change of coordinates are presented. 

Moreover, a family of Abel polynomial systems is also considered. 
By investigations until degree 10 we prove the existence of a unique 
C^ ' isochronous center. 

These results are established using a computer implementation 
based on Urabe theorem. \j 



>< 

c5 1 Introduction 



We consider the planar dynamical system, 



— = x = X{x, y), — =y = Y{x, y), (1) 
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where {x,y) belongs to an open connected subset U C M^, X,Y & C^(?7, R), 
and k > 1. An isolated singular point p e f/ of ([I]) is a center if and only if 
there exists a punctured neighborhood V C U oi p such that every orbit in 
V^ is a cycle surrounding p. 

The period annulus of p, denoted Vp is the largest connected neighborhood 
covered by cycles surrounding p. The period function T : Tp — > M associate 
to every point (x, y) G Tp the minimal period of the cycle 7(x,j/) containing 
{x,y). 

We say that a center p is isochronous if the period function is constant 
for all cycles contained in Tp. The simplest example is the linear center at 
the origin O = (0, 0) given by the system x = —y, y = x. 

For a cycle 7 G Fp we denote by C(7) C U the open subset bounded 
by 7. We say that the period function is strictly increasing (decreasing) iff 
T(7i) < T(72), (T(7i) > T(72)) for all 71 and 72 such that C(7i) C ^(72). 

An overview of J.Chavarriga and M.Sabatini [1] present the recent results 
concerning the problem of the isochronicity, see also IHEj. 

The main purpose of this paper is the study of the Lienard type equation 

x + f{x)x^ + g{x)=0 (2) 

with rational / and g, or equivalently the study of its associated two dimen- 
tional (planar) system 

i = y \ 

y = -gix) - fix)yj 

The Lienard type equation (l2|) appear for the first time in M.Sabatini 
paper |11) . when the sufficient conditions of the isochronicity of the origin O 
for the system (I3|) with / and g of classe C^ are given. 

In the analytic case, the necessary and sufficient conditions for isochronic- 
ity are given by A.R.Chouikha in [6], where the particular case of system 
system ^ 

x = -y + hx^y 1 

y = X + aix^ + a^y'^ + a^x^ + a^xy'^ J 

is studied. In this system as well as in all other considered systems, all 
parameters are real. All the values of the parameters for which the above 
system has the isochronous center at the origin are found. 



In [7] a similar result was obtained for more general system 

X = —y + axy + bx^y 1 

y = X + aix^ + asy^ + a^x^ + a^xy'^ J 

The aim of this paper is to extend investigations made in [6l [7] for sys- 
tems with higher order perturbations of the linear center x = —y, y = x. 
We investigate the practical applicability and the limitations of the method 
developed in the cited papers for more complicated systems. 

First let us consider the following particular case of ([3]) which is more 
general then ([5]) 

X = -y + bi^iyx + b2,iyx^ + 63,1^2;^ 1 , . 

, 2, 3, 2, 2, 22, 4f \^) 

y = X + a2flX + a^fiX + ao,2Z/ + 0'i,2xy + 02,2^; y + a^^x J 

Because of computational complexity, we select for investigation two sub- 
families (first one 61^1 = 03^0 = 0, second one bi^i = 62,1 = 0) of the above 
system which have the codimension two in the parameter space. 

In section 3, for the selected families we found all the parameters values 
for which the origin is an isochronous center. Thanks to this, among other, 
we found three new additional isochronous cases of linear center perturbed 
by homogeneous polynomial, which are not covered by the classification es- 
tablished by Chavarriga, Gine and Garcia in [2]. For these three isochronous 
centers we give the explicit form of the first integral and the linearizing change 
of coordinates. 

In the Section 4, an other particular case of the system ([3]) is considered, 
namely the fifth degree homogeneous polynomial perturbation of linear center 

„4 



X = —y + ayx 

y = X + bx^y^ + ex 



3„,2 , _5 ( (7) 



We found all the parameters values for which the center at the origin is 
isochronous (two families). The explicit form of the first integral and the 
linearizing change of coordinates are given for them. These systems are not 
contained in the Chavarriga et al. classification in [3]. 

In the last section, we investigate the following particular Abel polynomial 
system 

x = -y 

(8) 



y = J2Pk{x)y' 



^k{X)C 
fc=0 



where Pk{x) := a^x and a^ G M, for /c = 0, . . . , n. This Abel system is also a 
particular case of (I3|), and hence we can use the C-algorithm to investigate 
its isochronicity. 

Volokitin and Ivanov [12] proved that for n = 3 among systems of the 
form ([HD with arbitrary polynomials Pk{x) G W[x], there is only one family 
of isochronous centers. For Pk{x) = a^x, this family reduces to exactly one 
system. Namely, the following one 

i = —y \ 

y = x{l + yf] ^^^ 

In the cited paper Volokitin and Ivanov formulated the problem, which re- 
stricted to Abel equations of the form ([8]), can be stated as follows. Do exist 
among systems ^ isochronous ones with n > 4? We give a partial negative 
answer to this question showing that for 4 < n < 9 among systems ((HI) there 
is no an isochronous one. 



2 Efficient algorithm for computing 
necessary conditions of isochronicity 

2.1 About isochronous centers 

We collect now the results concerning Lienard type equation ((2|) (or its asso- 
ciated planar system (l3|)) which will be used later. 
Consider the Lienard type equation 

X + f{x)x^ + g{x) = 0, 

where / and g are C^ class functions defined in a neighborhood J of e M. 
Let us define the following functions 

PX PX 

F{x) := / fis)ds, 0(x) := / e^^'^ds. (10) 

Jo Jo 

When xg{x) > for x 7^ 0, define the function X by 

\x{xf = j\{s)e'^^^)ds. (11) 



Theorem 1 (H.Poincare). The planar system ^ with analytic data has an 
isochronous center at the origin if and only if for some analytic change of 
variables u = u{x, y) = x + . . . , v = v{x, y) = y + . . . the system ([T]) reduces 
to ii = —kv, V = ku, where fc G M, /c 7^ and . . . denotes the higher order 
terms. 

For more details see ITOl. 



Theorem 2 (Sabatini,[n]). Let f,ge C^{J,R). If xg{x) > for x ^ 0, 
then the system (JS]) has a center at the origin O. When /, g are analytic , 
this condition is also necessary. 

When f,g& C^{J,'R), the first integral of the system ([3]) is given by the 
formula 

/(x, x) = 2 r g{s)e^^^'Ms + (xe^^"))' (12) 

Jo 

Theorem 3 (Chomkha,[6]). Let f , g be functions analytic in a neighborhood 
J ofO, and xg{x) > for x ^ 0. Then system ([3l) has an isochronous center 
at O if and only if there exists an odd function h which satisfies the following 
conditions 

the function (f){x) satisfies 

j-Xix) 

0(x) = X{x) + / h{t)dt, (14) 

and X(x)0(x) > for x 7^ 0. 

In particular, when f and g are odd, then O is an isochronous center if 
and only if g{x) = e~^^^^(p{x), or equivalently h = 0. 

The function h is called Urabe function. As a corollaries of the above 
theorem one has 

Theorem 4 (Chouikha,[6j). Let f, g be functions analytic in a neighborhood 
0/0 G M, and xg{x) > forx ^ 0. If g'{x) + g{x)f{x) = 1 then the origin O 
is isochronous center of system (l3|) and its associated Urabe funtion h = 0. 

Theorem 5 (Chouikha,[6]). Let f , g be functions analytic in a neighborhood 
J ofO, Consider the system (l3|) having a center at the origin O. Let 

Sif.g) = 5/^(0) + 10/(0)/(0) + 8f (0) - 3/'(0) - 6r(0). 



Then the following holds: 

(a)- S{f,g) > then the period function T increases in a neighborhood ofO. 
(h)- S{f,g) < then the period function T decreases in a neighborhood ofO. 
(c)- If ([3]) has an isochronous center at then S{f,g) = 0. 

2.2 Chouikha algorithm 

The above Theorem [3] leads to an algorithm, first introduced by R. Chouikha 
in [6] (see also[7]), in what follow called C-algorithm, which allows to obtain 
necessary conditions for isochronicity of the center at the origin for equa- 
tion (ED- 

Below we recall basic steps of the algorithm. 

Let h be the function defined in the Theorem^ and u = 0(x). We assume 
that function (f) is invertible near the origin. 

where now X is considered as a function of u. Our further assumption is 
that functions f{x) and g{x) depend polynomially on certain numbers of 
parameters a := (ai, . . . , Op) e MP. 

By Theorem m if the system ^ has isochronous center at the origin, then 
the function h which is called the Urabe function, must be odd, so we have 

oo 

h{X) = J2c2k+iX"'+\ (16) 

A;=0 

and moreover, 

g{u) = g{x)e^^'^\ where x = (j)~^{u). (17) 

Hence, the right hand sides of (IT5|) and (ITTl) must be equal. Hence, we 
expand the both right hand sides into the Taylor series around the origin 
and equate the corresponding coefficients. To this end we need to calculate 
k-ih derivatives of flTSl) and flTTl) . 

For (ITSl) . by straightforward differentiation, we have 

d^~g{u) _ d [d^~^~g{u)\ dX 



du^ dX \ du^ ^ ) du 



Using induction, one can show that for ( TTTl) we obtain 



q^ = e(^-^-)^(^)^.(x), (19) 

where 5*^(0;) is a function of f{x),g{x) and their derivatives. 

Therefore to compute the first m conditions for isochronicity of system ([2]) 
we proceed as follows: 

1. we fix m and write 

m 
h{X) = Y, C2k-lX^'''' + 0(X2™), C := (Ci, C3, . . . , C2m-l), 
k=l 

2. next, we compute 

do 
'^k ■= Tt(0)' ^^^ = 'S'fc(O) 

for fc = 1, . . . , 2m + 1. Note that those quantities are polynomials in a 
and c. 

3. by the Theorem [3] we obtain equations v^ = Wk ioi k = 1, . . . , 2m + 1. 

It appears that we can always eliminate parameters c from these equa- 
tions. In effect we obtain a certain number of polynomial equations 
Si = S2 = S3 = . . . = sm = with p unknows Oj. These equations gives 
m necessary conditions for isochronicity of system ([2]). 

For more details see [H [7]. 

2.3 The choice of an appropriate Grobner basis 

Let us consider the following system 



X = -y + biAyx + . . . + 6„,_i,iyx" ^ 

y = x + a2flx'^ + ao,22/^ + . . . + a„-2,2a;''"^?/^ + a„,oa;" 

which is reducible to the equation ((21) with 



(20) 



ao,2 + 61, 1 + . . . + (an-2,2 + (n- l )bn-i,i)x'' ^ 
1 - 61,1a; - ... - 6„_i,ia;" 



/(^) = -. — J—: 1 z;^i ' (21) 
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g{x) = {x + a2flX^ + ... + a„,oa;")(l - &i,ia; - ... - 6„_i,ix" ^). (22) 

In this paper we have investigated the two types of high degree polynomial 
perturbations, homogeneous and non-homogeneous ones. It seems that C- 
algorithm is efficient for computing isochronicity necessary conditions for 
higher degree homogeneous perturbations. In this case system (l20l) reduces 
to the following one 



fc-i 



y = x + ak^2,2x''~^y^ + ak,ox'' J ^ ' 

where k E {2, . . . ,n}. 

For this homogeneous polynomial perturbation of the linear center, C- 
algorithm generate homogeneous polynomial equations in the parameters 
O'k-2,2,0'k,o and bk-i,i. Solving these polynomials, gives all the parameters 
values for which the real polynomial differential system ( !23ll is isochronous 
at the origin. 

However, for the non-homogeneous perturbation case, C-algorithm gener- 
ate non-homogeneous polynomial system. Moreover, non-homogeneous per- 
turbations depend on a bigger number of parameters. 

We note that for n = 3, C-algorithm succeeds to establish isochronicity 
criteria, however for n = 4 the obtained polynomials from the algorithm are 
much more involved. For example, for the system ( l20l ) with n = 4 reduces to 



X = -y + bi^iyx + b2,iyx'^ + 63,11/a;^ 

y = X + a2fix'^ + asfiX^ + 00,22/^ + ai,2xy'^ + a2,2x'^y'^ + a^flx'^ 

its associated first two non-zero polynomials obtained by applying C-algorithm 
are the followings 

P2 = 3 62,1-3 ai^2+bi/-a2fibi,i-9 ^3,0+4 ao,2^-5 &i,iao,2+10 a2/+l0 02,0^0,2, 

(24) 



Ps = 72 62,1^ + 396 02,0^1,101,2 + 90 61,100,201,2 + 36 61,102,2 + 324 63,100,2 

- 36 62,101,2 - 468 02,061,162,1 + 612 02,062,100,2 - 4116 61,102,0^00,2 

+ 10802,063,1 - 54003^062,1 - 324 04,061,1 + 1566o3_o6i,iOo,2 - 28802,002,2 

- 459 03,061,1^ - 1296 04,oOo,2 - 306 62,161,100,2 + 1428 02,061,1^00,2 

+ 15362,161,1' - 117 61,1 V,2 - 191 02,061,1^ + 180o2,oOo,20i,2 + 43 6i,/ 

- 2319 02,061,100,2^ - 289 61,1^00,2 - 360 00,202,2 - 3601,2' - 17162,100,2' 
+ 513 03,000,2' + 537 6i,i'oo,2^ + 351 Oo,2'oi,2 - 271 6i,iOo,2^ + 542 02,000,2^ 
+ 756 O2,o03,o«o,2 + 2268 02,003,061,1 - 20 00,2^ + 1120 02,0"^ + 798 6i,i'o2,o^ 

- 2240 61,102,0^ - 1512 02,004,0 + 100802,0^62,1 - 252o2,o'oi,2 

+ 1806o2,o'oo,2^ + 2240 02,0^00,2 

To solve the first nine non-zero obtained polynomials requests high perfor- 
mance computer and the standard accessible computer algebra systems for 
solving polynomial equations are not able to find a solution. 

For solving polynomial equations the Grobner bases are used. It is well 
known, see [I3|, that the form and the size of the Gobner basis of a polynomial 
ideal depends strongly on a choice of monomial ordering. A bad choice of 
the monomial ordering can be a main reason why the Grobner basis cannot 
be practically determined. 

Our basic observation concerning algebraic structure of polynomial equa- 
tions which give necessary conditions for the isochronicity in a case of non- 
homogeneous perturbations is following. Although the polynomials are not 
homogeneous, a careful analysis shows that they are quasi-homogeneous. In 
fact, one can notice that polynomial P2 given by (j24l) is homogeneous if we 
give weight 2 for 62,1, 01,2 and 03,0, and weight 1 for the remaining variables. 
More importantly we can find such a choice of weights for which all the 
obtained polynomials are homogeneous. 

The above observation shows that our main problem, i.e., finding a Grob- 
ner basis, concerns as a matter of fact, finding a Grobner basis of a homoge- 
neous ideal. It is well know, see [131 P- 466], that homogeneous Grobner bases 
have many 'nice' properties which make them extremely useful for solving 
large and computationally demanding problems. 

In fact, for non-homogeneous case of (1201) . the use of weighted degree 
gives a homogeneous Grobner base. 

To incorporate our observation into the C-algorithm we choose a new 
parametrization for the problem. First, we observe that all polynomials which 



are obtained by means of the C-algorithm are homogeneous if we choose the 
following weights 

1. z + j — 1 for parameters a^j and bij 

2. 2z + 1 for C2i+i. 

Knowing this we introduced new parameters Aij, Bij, and C2i+i putting 

^:y"' = «M Bl^'~' = k^^ C^l = c,.^, (25) 

After this reparametrization system (|20|) reads 

x = -y + Bi^iyx + ... + BUzl.yx''-' 

y = x + A2flX^ + Ao,2Z/' + . . . + AlzUx^'-'^y' + AlZ^'x 



2 , A „.2 , , AU-l ri-2 2 , An-l^n ( \^^l 



As in the case of isochronous center the Urabe function is odd, we search it 
under the form 

oo 

h{X) = J2 Cgfc+i^^^^^ = ^1^ + ^1^^ + ^5^^ + ^7^^ + • • • (27) 

fc=0 

We emphasize that from the isochronicity conditions for fl26|) . expressed in 
terms of its parameters, it is easy to reconstruct the parameters values for 
which the system f l20l) admits isochronous center at the origin, by a simply 
use of (ESD. 

Fact, that the described reparametrization gives rise homogeneous equa- 
tions, allows to reduce the number of the parameters appearing in (l26l) by 
one. First, we assume A2.0 = 0, and then solve the isochronicity problem for 
system (l26l) under this assumption. Next, for ^2,0 7^ 0, we apply on ( !26l ) the 
following change of coordinates 



{x,y)^--{x,y) (28) 

^2,0 



We obtain 



^ =-^+(ft^^^ + - + te' ^^""' 



n-l 



n—1 / , \ n—1 



^ + ^^+(tfh/' + - + Hrf ^"-V+ tf) ^" 
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(29) 



Hence, without loss of generality we can put ^2,0 = 1, and find the parame- 
ters values for which the center is isochronous. 

We note that for an arbitrary k E N, the problem of the isochonicity of 
the center for homogeneous perturbations of the form ([23l) reduces to solve 
a number of polynomial equations in 3 parameters. 

Recall that linear center perturbed by homogeneous polynomial, was in- 
vestigated by W.S. Loud in [8] for the quadratic case, and in [6] the author 
find Loud results by the described algorithm, see also [4]. 

Homogeneous perturbations was also studied by Chavarriga and cowork- 
ers. For the fourth and fifth degree homogeneous perturbations, see [21 [3], 
where the homogeneous perturbations different from those studied in the 
present paper are considered. 

3 Fourth degree perturbations 

Let us consider the following system 

X = -y + bi^iyx + 62,12/2:^ + bs^iyx^ 1 .^^. 

y = X + 02,0^^ + 03,0^^ + ao,22/^ + ai,2xy'^ + 0.2,23^^2/^ + a4.,ox* J 

For the system (1301 ) having a center at the origin O, the following lemma 
gives a montonicity criteria for its period function 

Lemma 1. Let 

S = lOalo + 10a2,oao,2 — 3ai,2 + 362,1 — &i, 1^2,0 + ^1,1 — 803,0 + 4ao,2 ~ 56i,iao,2 

If S > (S < 0) then the period function of ( l30l) is increasing (decreasing) 
atO. 

Proof. System ( l30l ) reduces to the Lienard type equation ([2]), with 

„. X _ Qo,2 + fei,i + (ai,2 + 262,1)3^ + (^2,2 + 3&3,i)x^ 1 

^^'''~ 1 - 6i,ix - 62,1x2 - 63,1x3 [ (31) 

g{x) = (x + a2fix'^ + a3fiX^ + a4,oa;^)(l - 61,1a; - 62,12;^ - 63,ix^) J 

then, we establish only three iterations (derivatives) of the C-algorithm 
given in Section 2, after elimination of Ci of the Urabe function from the 
second derivatives, substitution in the third polynomial gives 5". Using the 
Theorem [HI (See Corollary 2.8 of [6]), we prove the result. D 

11 



Because of computational complexity, we select for isochronicity investi- 
gation, two sub- families of the system (1301 ). which have the codimension two 
in the parameters space. 

3.1 First family 

Let us assume 61,1 = 03,0 = 0, in this case (1301) reduces to the system 

X = -y + h2,ix^y + h^iyx'^ 1 ,^^. 

y = x + a2flx'^ + ao,2y^ + ai,2xy'^ + «2,2a;^?/^ + 04,03^^ J 

For this system having a center at the origin O, we give isochronicity 
necessary and sufficient conditions depending only on the following seven 
real parameters 62,1, ^3,1, ^2,0, 00,2, «i,2, ^2,2, 04,0- 

Theorem 6. The system (l32l) has an isochronous center at O if and only if 
its parameters satisfy one of the folowing conditions 

1. 02,0 = O0,2 = ^2,1 = Cil,2 = 0, &3,1 = -4a4,o/3, 02,2 = -16(3-4,0/3 

2. 02,0 = ao,2 = ^2,1 = ai,2 = 0, 04,0 = -&3,i/2, 02,2 = ^3,i/2 

3- 02,0 = '^4,0 = (^0,2 = ^2,1 = fll,2 = 0, 02,2 = &3,1 

which give the homogeneous perturbations and remaining ones which give the 
non-homogeneous : 

4- 0^2,0 = 0^4,0 = '^0,2 = 0, 02,2 = ^3,1) ^2,1 = ^1,2 

5- Oo,2 = — 2o2,o, 0-1,2 = 8*22 0/3, 02,2 = ~8o2,o/3) 

«4,o = 0, 63,1 = -4o|o/3, &2,i = 2o|o/3 

6. oo,2 = -2o2,o, ai,2 = 8020/3, 02,2 = -16o2,o/21, 
04,0 = -4o3 0/21, 53,1 = -4o|o/21, &2,i = 2o2 0/3 

'^- Ct2,0 T" 0? '^0,2 = — 202,0, Ct4,0 = 0, 02,2 = 263,1, 

ai,2 = ai,o(4 + &3,i/a^,o). hi = 4o(2 + ^3,1/^2,0) 
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In our investigations we have used Maple in its version 10. To compute 
the Grobner basis of the obtained polynomial equations in the ring of charac- 
teristic 0, we have used Salsa Software more precisely the FGb algorithm see 
[T4] . In this proof, we do not present the algorithm generated polynomials 
which are too long. 

Proof. We consider separately the two cases : homogeneous and non-homogeneous 
perturbations. 

Homogeneous perturbations 



In the homogeneous case system ( 1321 ) reduces to 



y = X + a2^2X^y^ + o,4:,ox^ 



which is reducible to (Ej) such that 

(363,1 + a2,2)x 



2 



fix) 



l-&3,ix3 ' } (34) 

g[x) = (x - h^ix^) (1 + a4_oX^) , 

C-algorithm gives three homogeneous perturbations of linear center, can- 
didate to be isochronous. 19 derivatives was essential to obtain the necessary 
conditions of isochronicity. We give explicitly the Urabe function h associ- 
ated to each system given by the algorithm. 

Case 1: the system ( !33l ) becomes 

" = -' -J'^'^'li ^ 4 ] (35) 

y = X - -faifiX^y^ + a^fiX^ J 

We can easily check that f{x)g{x) + g'{x) = 1, hence, following Corollary 2-7 
of [6], the system (i35l ) is isochronous and h{X) = 0. 

Case 2: the system (l33l) becomes 

X = -y + h^iyx 

. ^ ^3,1 2 2 ^3,1 4 } (36) 

y = X -\ X y X 

^ 2 2 



3 
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The computations yield to the coefficients of the Urabe function : 

ci = 0, C3 = -^, C5 = 0, C7 = 0, 

^i,i _. _. _3&i,i 

C9 — "TTT' "^11 ~ "^' "^la — u, Ci5 — , 

lo 256 

Let Mo = C3, Ml = cg, M2 = C15, M3 = C21 . . . We observe that for A; = 1, 2, 3, 4 

Uk+i _ &i.i(fe + i)(fe + i) 

Mfc 4(A; + 1)2 

It is then natural to conjecture that this is always the case. By series sum- 
mation we found the odd function 



oo 



hiX) = ^M^= = J2u,X' (37) 

U + h/x^ fc=i 



By direct computations we verify that the above h satisfy the equation (fT3|) . 
Thus we conclude that the obtained h is the Urabe function. 
This case is similar to the one found in Theorem 3 of [7]. 



Case 3: the system ( 1331 ) becomes 

X = -y + a2,2yx^ 
y = x + a2,2X^y'^ 



(38) 



also this equation has an isochronous center at O since f{x)g{x) + g'{x) = 1 
and h{X) = 0. 

We note that the case 3 is the case 4 with ai 2 = 0. 

Non-homogeneous perturbations 

For the non-homogeneous perturbation of the linear center, we add to 
C-algorithm the two tricks expanded in the last Section ( Homogenization 
and reduction of the dimension of the parameters space by 1 ). 
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We obtain, the last four theorem cases 4, 5, 6 and 7, when restricted to 
ao,2 = 1, satisfy assumptions of Theorem [H we show that each of the last 
four systems restricted to a2,o = 1 admit isochronous center at the origin and 
h = 0. Finally, by the rescaling 

{x,y) ^ ao,2ix,y) (39) 

we generalize the isochronicity for any 02,0 7^ since by change of coordi- 
nate (i39l ) the isochronicity is not lost. 

n 

To complete the analysis, to each isochronous center obtained by homo- 
geneous perturnation, using Theorem [2] we write explicitely the first integrals 
and thanks to the method described in ^0\ we compute the linearizing change 
of coordinates. 

A first integral of the system (l35ll is 

_ {alox^ + 2a4,o x^ + x"^ + y"^) ^ 
729 (3 + Aa^fiX'^) 
The following analytic change of coordinates 

X (1 + 04^02^^) 
"~ 3 (3 + 4 04,0^3)^/^ 

^ y_ 

3 (3 + 4 04,0 x3)^/^ 

transforms the system (l35ll into the linear system : 

ii = —V 
V = u 

A first integral of system ( l36ll is 

^® - 63,1x3 - 1 ' 

then we give the linearizing change of coordinates 

X 

u 



^-l+&3,ia;3 

y 
^-1 + 63,1x3 
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A first integral of (i38ll is computed 



Gl - T~i~< 3^2 



A linearizing change of coordinates is 

X 

u 



\/-l + 0,2,2 X^ 

y 

^-1 + a2,2X^ 

For the remaining cases 3,4,5 and 6 obtained by non-homogeneous pertr- 
bations, the first integrals and linearizing transformation can be obtained by 
the same way, but the computations are cumbersome and we dont report 
them there. 



3.2 Second family 

Consider system (l30l) . with h^i = 62,1 = 0. We obtain the seven parameter 
real system of degree 4. 

x = -y + h,wx' 1 

I 2, 2, 3, 2, 22, 4f \^^J 

y = X + a2flX + ao,22/ + 03,00; + ai^2xy + a2,2X y + 04,02; J 

For this system having a center at the origin O, we give isochronicity neces- 
sary and sufficient conditions depending only on these seven real parameters 

03,05 ^3,l5 02,0, Oo,2; Oi,2, 02,2; 04,0. 

System (l40ll reduces to the equation ([2]) with 

ao,2 + oi,2X + (a2,2 + ^h,i)x^ 
•/ ^ = ^ 1 3 ' 

g{x) = (x + O2,ox^ + O3,ox^ + a4,ox'^)(l - 63,ix^) , 



Theorem 7. The system ( l40l) has an isochronous center at O if and only if 
its parameters satisfy one of the folowing conditions 

1. 02,0 = oo,2 = &2,i = oi,2 = 0, 63,1 = -4a4,o/3, 02,2 = -16a4,o/3 
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2. 02,0 = ao,2 = &2,i = ai,2 = 0, 04,0 = -&3,i/2, 02,2 = &3,i/2 

3- 02,0 = '^4,0 = «0,2 = ^2,1 = Ol,2 = 0, 02,2 = ^3,1 

which give the homogeneous perturbations and remaining ones which give the 
non-homogeneous 

4- «2,o = -«o,2/2, 03,0 = «o,2(9 - V33)/48, ai,2 = ao,2(-l + y33)/16 

04,0 = 0, 02,2 = ao,2(-21 + 5y33)/64, 63,1 = a|]^2(-21 + 5 V33)/192 

5. 02,0 = -00,2/2, 03,0 = 00,2^(9 + V33)/A8, 01,2 = -0^,2(1 + V33)/16 
04^0 = 0, 02,2 = -00,2(21 + 5^)/64, 63,1 = -olj 2(21 + 5 V33)/192 

6- 02,0 = —00,2/2, 03,0 = 04,0 = 0, ai,2 = 00,2/2 
02,2 = oL/2, 63,1 = ag 2/4 



7. 02,0 ^ 0, ao,2 = -2a2,o, J^ = ^ - Ss (22868 + 468 v/3297)''^ 
+ ^ Vsm (22868 + 468 ^3297)^^^ - j^ \/22868 + 468 ^3297, 

12,2 _ 344 _ 4720 V3297 i 96896 1 



-'2,0 



3549 273 /'oocco , ^co /^t^T^^^/s i 3549 3 



(22868+468 V3297) '^^^^ ^22868+468 v^3297 



+§- 3^ v^^ - 2^^ (22868 + 468 ^3297) "'^' , 

'^'-^ V 22868+468 v^3297 '^°'*^ ^ ^ 

_ 2 /-584+(22868+468 V3297)^''^+14 v^22868+468 ^3297 , 322 

Ol,2 = O2o(, 3/ jj 03,1 = -^, 

39 V 22868+468 \/3297 * 

a4,o _ 2150 _|_ 11926 1 2_ v^3297 

«2,o 10647 10647 ^22868+468 V3297 ^^ ^/22868+468 V3297 

1085248 (noQdQ 1 /1«C , /QOn7\^2/^ 160 v^3297 



^7^j^2^ (22868 + 468 V3297) ' - -^ -^73 

10647 V V ; 91 (^22868+468 ^329?)'^' 

Proof. The same method used in the last proof, is employed. We establish 
the seven cases given in the theorem. 

The three first cases belong to the family of linear centers perturbed by 
fourth degree homogeneous polynomial. Those are analyzed in the proof 
concerned by the first family . 

For each of the cases 4, 5 and 6 of the theorem restricted to ao,2 = 1, it 
is easy to verify that conditions of Theorem [4] are satisfied, also for the case 
7 restricted to 02,0 = 1 we check easily f{x)g{x) + g'{x) = 1. Then, by an 
appropriate change of coordinate, we show that the system ( !40l ) satisfying 
one of the cases 4, 5, 6 and 7 have isochronous center at the origin. D 
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Note that the polynomials issued from the 19 derivations and associated 
eliminations for the system (i30l) (with 9 parameters), exceed the authorized 
memory of ordinary computers (2 Go of Random Access Memory ) in com- 
putations of the Grobner basis by the known efficient algorithm FGb [T4] . 

4 Fifth degree homogeneous perturbations 

Let us consider the following system of degree 5. 

X = — V + ayx'^ I , , 

y = x + bx^y + cx^ \ 

The systems given in the following theorem are additional isochronous cases 
to those established by Chavarriga et al in [3]. 

Theorem 8. System (|4T]1 has an isochronous center at the origin if and only 
if it reduces to one of the following systems 



X = —y + ayx'^ 



3.2 4 .^ (42) 



y = X + 5ax y ax 

5 

x = -y + ayx'^\ 

3 2 } (43) 

y = X + ax y J 

Proof. We perform the C-algorithm with the functions / and g such that 

_ {b + Aa)x^ 
^^''^~ (l-ax^) 
g{x) = (x + cx^)(l — ax^) 

then we obtain isochronicity necessary conditions for the system ([4Tll.For the 
two systems ([42l) and (|43|) . it is easy to check that assymptions of Theorem [4] 
are satisfied then systems (l42l) and (j43|) are isochronous and their associated 
Urabe function /i = 0. D 

As before, for each obtained isochronous center, we write explicitely the 
first integrals and the linearizing change of coordinates. 
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For the system ( l42i ) a first integral is 

({-5 + 4:ax'^fx'^ + 25yA 
625 (—1 + ax^j 

We give the linearizing change of coordinates for the system (l42l) 

(—5 + 4 ax'') X 



u 



5 (-1 + ax^f'^ 



{-l + ax^f^ 
For the system fl43l) a first integral is 



^63} - — 



2 I 2\2 



ax^ 



We give the linearizing change of coordinates associated to ( 1431 ) 

X 



M 



v^-l + ax^ 

y 

v^-l+ax^ 



5 The period function for an Abel polynomial 
system 

This section is concerned with the following Abel system 

x = -y 

n 

y = ^Pk{x)y'' 

fc=0 

where Pfc(x) := a^x, a^ := 1, and a^ G M, for /c = 1, . . . , n. 



(44) 
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5.1 Reduction to Lienard type system 

The system (144|) can be written 

x = -y 

y = x{l + P{y)), 

with P{y) = aiy + 02?/^ + a^y^ + .... + a,,?/". 
Let define the functions X and as follow 

1 r s r ds 

-X(x)^ =: / — ds and 0(x) 



(45) 



2 ^ ^ 7o (l + i^(^)) ^' ' 'Jo a + Pis)) 

We show the following. 



Theorem 9. The origin O is a center for ( 1451) . 

Moreover, if a\ — 3a2 > {a\ — 802 < 0) then the period function of (l45!) 
«s increasing (decreasing ) at O. 

This center at O, is isochronous if and only if there exists an odd function 

h satisfying 

X 

1 + h{X) ~ "" 
and 



0(x) = X{x) + / h(t)dt 
Jo 

such that X(j){x) > for x 7^ 0. 



In particular, when P is an even polynomial then the origin is isochronous 
center if and only if P = Oa- 

Proof. From the symmetry criteria, the origin is a center for fl45ll . see [9l 
chap. 4]. We see that when -P(O) = 0, there exists an open connected interval 
Ji containing where 1 + P{y) 7^ 0. Then we can consider system in a 
neighborhood U of the origin where U = Ji x J2 with J2 a suitable open 
interval containing 0. So, making the following transformation : 

^ = ^ I (46) 

z = x{l + P{y)) ^ ^ 



^The last statement of the above theorem was communicated to the autor by 
A.R.Chouikha. 
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we obtain, 

y = z 



z = -y{l + P{y)) + z 



2 P\y) 



In effect, after renaming y as x and z as y we obtain the ((2l). The origin O 
is a center for ([2]) with 



(l + Pix)) 
g{x) =x{l + P{x)), 

For a fixed n > 2, we establish three iterations of the C-algorithm given 
in Section 2, after elimination of Ci of the Urabe function from the second 
derivation, substitution in the third polynomial gives a\ — 3a2. By Using the 
Theorem [5] we prove the monotonicity result. 
Let 



F{x)= f{s)ds = -\n{l + P{x)), 
Jo 

0(a;) = u= e^Wds = ' 



(1 + P(s)) 
Then we obtain 

g{u) = ^(x)e^(") = x(l + P(x))e"'"(^+^(")) = x 

Following Theorem [3l g satisfies 

, . X 



9[u) 



l + h{X) 



where u = X + H{X). 

For the particular case when P is even, it is easy to see that / and g are 
odd. We use the Theorem [HI that leads us to obtain x = X and P = 0. D 

The following paragraph is devoted to illustrate the last theorem by ex- 
ample. 



21 



5.2 Application to Volokitin and Ivanov system 

This section concerns the Abel polynomial system 

x = ~y 1 

i=0 ) 



with Pk eR[x], 0<k<n. 

For such system the origin O is not always a center. 

Theorem 10 (Volokitin and Ivanov, [E])- The polynomial Abel system 

x = -y 

y = {x + a^x^){l + h{x)yf + 3axy{l + h{x)yY - h'{x)y^ 



(48) 



with an arbitrary number a G M and an arbitrary polynomial h G ]R[x] has 
an isochronous center at O. 



Let us consider the Abel system ( 1441 ) with n = 9 : 

x = -y 

9 

y = x + '^ Oixy' 



(49) 



i=l 



with afc G M, 1 < A; < 9. As follows from the symmetry criteria, the origin O 
is always a center for (l49l) . see [9l chap. 4]. 



Theorem 11. On/y m t/ze case 

x = -y 



a\ 2 (^ 



L..3 (50) 



y = x + aixy + — xy + — xy 

the system (l49l l /zas an isochronous center at the origin O . That is, only in 
the case when the system ( l49l) belongs to the Volokitin-Ivanov class (i48l) . 



Proof. We perform C-algorithm with 

„. . P'(x) 

fix) = 



(l + Pix)) 
g{x) = x{l + P{x)) 
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where P{x) = '^i^iCLiX\ We obtain the unique one-parameter family (ISOll . 
and computations gives the Urabe function 

ciiX kiX 

h[X) = 



3 r (51) 



3 ^fWTk^ 

with ki = — ai/3, k2 = 1, k^ = 0. 

It is interesting to note that this Urabe function is of the same nature 
that the Urabe function ([371) of the system (l36ll . (cf. also the proof of the 
Lemma 3.4 of [6| as well as Section 3 and 4 of [7].) D 

We transform system (iSOl ). by the following change of coordinates 

^ ~ 3 

^ ~ 3 ■ 

Then, after renaming ^ as a; and C as y we obtain the Volokitin and Ivanov 
system [12], p. 24], which is a particular case of (1481) 

x = -y 

y = xil + yY. 

In [I2], they prove that O is an isochronous center of (jSTl ) showing that it 
commutes with some transversal polynomial system, but dont provide its 
first integral. 

Thanks to our rescaling (146]) . we determine a first integral of (ISTl) 

Indeed, let P{x) = 3x + 3x'^ + x^. The system ( ISTl ) reduces to the system ((3|) 
with 

^' ' (i + pw) 

g{x) =x(l + P(x)). 
By Theorem [2] we have 

/(x,x) = 2 /%(s)e2^(^)rfs + (xe^("))2. 
23 



is a first integral of ([3]) which, in our case, reduces to 



By the reciprocal of the rescaling fl46ll . we obtain the first integral of the 
system (ISTll 

Unfortunately we are unable to find explicitely the linerizing transforma- 
tion for the system (ISB . 

In the light of Theorem [H], it is natural to ask if the system ( l50l ) is the 
unique system with isochronous center at the origin inside the family ( 1441) . 

Even for the system (j44l) whith n = 10, our actual computer possibilies 
are not sufficient to give an answer. 
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